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Abstract 

In any reaction-diffusion system of predator-prey models, the population densities of 
species are determined by the interactions between them, together with the influences 
from the spatial environments surrounding them. Generally, the prey species would die 
out when their birth rate is too low, the habitat size is too small, the predator grows too 
fast, or the predation pressure is too high. To save the endangered prey species, some 
human interference is useful, such as creating a protection zone where the prey could cross 
the boundary freely but the predator is prohibited from entering. This paper studies the 
existence of positive steady states to a predator-prey model with reaction-diffusion terms, 
Beddington-DeAngelis type functional response and non-flux boundary conditions. It is 
shown that there is a threshold value do which characterizes the refuge ability of prey 
such that the positivity of prey population can be ensured if either the prey’s birth rate 
satisfies 0 > 6q (no matter how large the predator’s growth rate is) or the predator’s 
growth rate satisfies n < 0, while a protection zone Do is necessary for such positive 
solutions if d < do with ^ > 0 properly large. The more interesting finding is that there 
is another threshold value 9* = d*(^. Do) < do, such that the positive solutions do exist 
for all d S (d*,do). Letting fi —>■ oo, we get the third threshold value di = di(Do) such 
that if d > di(Do), prey species could survive no matter how large the predator’s growth 
rate is. In addition, we get the fourth threshold value d* for negative fj, such that the 
system admits positive steady states if and only if d > d*. All these results match well 
with the mechanistic derivation for the B-D type functional response recently given by 
Geritz and Gyllenberg [A mechanistic derivation of the DeAngelis-Beddington functional 
response, J. Theoret. Biol. 314 (2012) 106-108]. Finally, we obtain the uniqueness of 
positive steady states for p properly large, as well as the asymptotic behavior of the 
unique positive steady state as ^ oo. 
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1 Introduction 


Biological resources are renewable, but many have been exploited unreasonably. Nowa¬ 
days, some species cannot survive in their habitat withont human intervention. Snch inter¬ 
ventions have included establishing banned fishing areas and fishing periods to cope with over¬ 
fishing in fishery production, and setting up nature reserves to protect endangered species. 
These phenomena are usually described via diffusive predator-prey models, where the popu¬ 
lation evolution of the species relies on the interactions between predator and prey, as well 
as the influences from the spatial environments surrounding them. Naturally, prey species 
would die out when the prey’s birth rate is too low, the habitat size is too small, the preda¬ 
tor’s growth rate is too fast, or the predation rate is too high. To save the endangered prey 
species, various human interferences are proposed such as creating a protection zone where 
the prey could cross the boundary freely but the predator is prohibited from entering. Refer 
to the works on protection zones by Du et al for the Lotka-Voltera type competition sys¬ 
tem [7], Boiling n type predator-prey system [9], Leslie type predator-prey system [8], as 
well as predator-prey systems with protection coefficients [10]. Oeda studied the effects of a 
cross-diffusive Lotka-Voltera type predator-prey system with a protection zone [22] . A cross- 
diffusive Lotka-Voltera type competition system with a protection zone was investigated by 
Wang and Li [27]. Zou and Wang studied an ODE model of protection zones, where the sizes 
of the protection zones are reflected by restricting the functionals’ coefficient for the predator 
[31]. Recently, Cui, Shi and Wu observed the strong Allee effect in a diffusive predator-prey 
system with protection zones [4]. 

In this paper, we study the steady states to the following diffusive predator-prey system 
with Beddington-DeAngelis type functional response 


ut — diAu = u{6 — u — 

Vt — d2Av = v{fi — V + 
du 

dn ’ 
dv 

dn 

u{x,0) = uo{x) > (^)O, 
^ n(x,0) = vo{x) > (^)O, 


a{x)v 


1 -I- mu + kv^’ 
cu 

1 -|- mu + kv ' 


X G D, 

t > 0, 

X G Di, 

t > 0, 

X G cID, 

t > 0, 

X G cIDi 

, t > 0 

X G D, 


X G Di, 



( 1 . 1 ) 


where D is a bounded domain in {N < 3) with smooth bonndary dQ, Dq d D with 
di^o smooth, = D\Do, constants di,d 2 , 0 ,c,m,k > 0, // G M, ^ is the outward normal 
derivative on the boundary, and 


o(x) 


0, x € Do, 

a, X G Di. 


( 1 . 2 ) 


The fact of a(x) = 0 in Dq implies that no predation could take place there. 


2 






Eq. (1.1) is a reaction-diffusion system of species u and u, and the dynamical behavior of 
species would be determined not only by the mechanism of the functional response between 
u and u, but also by the interaction between their reaction and diffusion. Here prey u and 
predator v disperse at rates di and d 2 , and grow at rates 9 and //, respectively. The prey 
is consumed with the functional response of Beddington-DeAngelis type in and 

contributes to the predator with growth rate in Hi. Non-flux boundary conditions 

mean that the habitat of the two species is closed. The B-D type functional response was 
introduced by Beddington [1] and DeAngelis [5]. Refer to [1, 5, 6] for the background of the 
original predator-prey model with B-D type functional response. Guo and Wu studied the 
existence, multiplicity, uniqueness and stability of the positive solutions under homogeneous 
Dirichlet boundary conditions in [13], as well as the effect of large k in [14]. Chen and Wang 
established the existence of nonconstant positive steady-states under Neumann boundary 
conditions [2, 23]. 

In particular, a mechanistic derivation for the B-D type functional response has been 
given by Geritz and Gyllenberg in [12] recently, where predators v were divided into searchers 
vs with attack rate a and handlers vh with handling time h, while preys u were structured 
into two classes: active preys up and those prey individuals up who have found a refuge 
with total refuge number b and sojourn time r. In these terms, the parameters in B-D type 
functional response of (1.1) can be understood as that m = ah reflects the handling time of 
Vh, and k = hr describes the refuge ability of the prey. 

The prey’s refuge may come from its aggregation, reduction of its activity, or places 
where its predation risk is somehow reduced [26]. Dynamic consequences of prey refuges 
were observed by Gonzalez-Olivares and Ramos-jiliberto with more prey, fewer predators 
and enhanced stability [11]. On the other hand, refuges from species usually cost the prey in 
terms of reduced feeding or mating opportunities [26], and hence their population could not 
be very large. In contrast, the protection zones, as refuges from humans, always benefit the 
endangered species. Refer to [15, 16, 21, 25, 28, 29, 30] for more backgrounds on prey refuges 
and their affections. In this paper, we will show the effect of the prey’s refuge and the size 
of the protection zone on the coexistence and stability of the predator-prey system with B-D 
type functional response. The results obtained here observe the general law that refuges and 
protection zones benefit the coexistence of species [11, 26, 31]. 

Since the model (1.1) contains different coefficients a{x) and c in the B-D type functional 
response terms for u and v respectively, without loss of generality, suppose di = ^2 = 1 for 
simplicity. The steady-state problem corresponding to (1.1) takes the form 


—Au = u{9 — u — 

—Av = v{fi — V + 
du 


a{x)v 


dn 


dv 

do. on 


1 -|- mu + kv' 
cu ^ 

1 -|- mu + kv' 

= 0 . 


in H, 
in Hi, 


(1.3) 


an. 


Denote by Xi{q) the first eigenvalue of —A -|- q over H under homogeneous Neumann 
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boundary conditions with q = q[x) G L°°{Vt). The following properties of Ai(g) are well 
known: 


(i) Ai(0) = 0; 

(ii) Ai(gi) > Xi{q2) if qi > q2 and qi ^ q2] 

(iii) Xi{q) is continuous with respect to g G L°°{Q). 
Define 


Do) = Xi{q{x)), ^0 = ^, 6'i(Do) = Xi{qo{x)), 


(1.4) 


with 


a(x)/r 


0, X G Do, 
6*0, X G Di. 


(1.5) 


Denote by Uq the solution of the scalar problem 


Ou 

—Au = u{9 — u — qo(x)) in Q, — = 0 on Oft. 

on 


( 1 . 6 ) 


Tq Tq (jP'dx 

Due to 01 = inf - ^ , 2 , —^-, the properties (i)~(iii) of Xi(q) imply 

the following lemma: 


Lemma 1.1 0*(/x,Do) is strictly increasing with respect to fj, and decreasing when Do enlarg¬ 
ing, 0*(O,Do) = 0, 0*(/x,Do) < 00, lim/,^oo0*(M,^o) = 0i(Do) < lim|Qj|^o ^ 0 ) = 
0, lim|iig|^o^'*(/^,^^o) = T^- □ 

Biologically, we are interested in the positivity of the prey u in the diffusive predator-prey 
model (1.3). We state the main results of the paper one by one as follows. 

Obviously, either large 0 or small /r benefits the prey u. In the first theorem, we give 
two sufficient conditions for keeping the prey positive without protection zones. 


Theorem 1 If 0 > 6 q or // < 0, then the positivity of u would be ensured automatically 
without any protections zones. 

The next theorem implies that a suitable protection zone guarantees the existence of 
positive solutions to (1.3) under 0 < 0o with p > 0. 

Theorem 2 Suppose p > 0. If 0*(^, Dq) < 0 < 0o, then Eq. (1.3) has at least one positive 
solution. Furthermore, if 9 < 0*(//, Dq) with m < , then Eq. (1.3) has no positive 

solutions. 


In the third theorem, we give a necessary and sufficient condition for the coexistence of 
u and V under p G (— ;^,0]. 
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Theorem 3 Suppose — ^ < /i < 0. Then Eq. (1.3) has at least one positive solution if and 
only if 6 >e, = > 0. 

Remark 1 Since ^o) = 0 by Lemma 1.1, for any 6 > 0 and p. > 0 fixed, 

the key condition 9 > in Theorem 2 can be realized by enlarging the size of the 

protection zone Qq = So does the condition 6 > 6i m the following Theorem 4. 

Remark 2 Theorem 1 shows that no protection zones are necessary for the positivity of u if 
p < 0. It is known by Theorem 3 that in addition to the positivity of u, the positivity of v can 
be ensured also if the death rate of the predator v is not too high with p G (—^,0] C (—oo,0] 
and the birth rate of the prey u is properly large such that 9 > 6^,. 

Finally, the last theorem says the positive solutions of (1.3) are in fact unique if 9 is even 
larger than 9i under large p, and determines the asymptotic behavior of the unique positive 
solution as /i —>■ oo. In fact, from Lemma 1.1 and Theorem 2 that if 0 > 9i, prey species 
could be alive no matter how large the predator’s growth rate is. 

Theorem 4 If 9 > 0i(no)? then there exists p* > 0 such that the positive solution of (1.3) 
is unique and linearly stable when p > p*. Furthermore, the unique positive solution satisfies 
{u,v — p) ^ {Ug^qQ,0) uniformly on and fli, respectively, as p ^ oo. 

This paper is arranged as follows. In the next two sections, we give the proofs of The¬ 
orems 1-3 and Theorem 4, respectively. The last section is devoted to a discussion of the 
obtained results, by analyzing them with the mechanistic derivation for the B-D type func¬ 
tional response in [12]. 


2 Existence of positive solutions 

At first we deal with the proof of Theorem 1. 

Proof of Theorem 1. Assume p < Integrate the second equation of (1.3) over fli to 
get 


0 = / v[p — V + 

J 


cu 


1 -|- mu + kv 


dx, 


and hence 


0 < / v'^dx = / v{p 




cu 


1 -|- mu + kv 


)dx < {p + 


-) j 


vdx < 0. 


This concludes u = 0, and so u satisfies 


du 


—Au = u(9 — u) in fl, = 0 on dH. 

on 


( 2 . 1 ) 


Obviously, (2.1) admits the solution u = 0 > 0. 
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The desired result for < /i < 0 is substantially concluded from Theorem 3. Indeed, 
the subcase of 0 > 0* is covered by Theorem 3, while for 0 < 0*, it can be found in the proof 
of Theorem 3 that u = 0, and so u = 0 > 0. 

Next consider the first equation of (1.3) with 6 > 9q. It is easy to know that < 

< f = ^0 for u > 0. Thus, for any v{x) > 0, there is 9q € (0, 9q) such that 

dix^v ~ du 

—Au = u{9 — u — - - - —) > u{9 — 9o — u) in Q, — = 0 on dfl. 

1 + mu + kv on 

This ensures that u > 9 — 9q > ^. □ 

We need some preliminaries represented as lemmas and propositions for the proof of 
Theorem 2, and begin with two known results on the maximum principle and the Harnack 
inequality. 

Lemma 2.1 (Maximum Principle [19]) Let g € C{Ll x M), w € (17^(12) f) (7^(12), where 12 is 
a bounded domain in with smooth boundary. 

(a) If Aw + g{x, w) < 0 in Ll, ^ > 0 on (912 and minjjrc = w{xo), then g{xo,w{xo)) < 0. 

(b) If Aw + g{x,w) >0 in Ll, ^ <0 on dLl and max^u) = w(xo), then g{xQ,w{xQ)) > 0. 

□ 


Lemma 2.2 (Harnack Inequality [20, 17]) Let f G L^iLl) with p > max{Y,l}, and w be 
a non-negative solution of Aw + f{x)w = 0 in a bounded domain 12 C with smooth 
boundary under homogeneous Neumann boundary condition. Then there exists a positive 
constant C = C{p,N,Ll, ||/||lp(o)) such that 

raayiw < Cramw. □ 
o f! 


The following a priori estimates are easy to get. 

Lemma 2.3 Let {u,v) be a nontrivial non-negative solution of (1-3). Then 


D < u < 9, < u < + 


c9 


l + m9 + k^i+^ ll«llci.-(o) + Ibllci-(Oi) < 
with = max{/z, 0}, a G (0,1), C = C{9, p, LIq) > 0. 

Proof. Suppose u{xq) = max^u(x) > 0. By Lemma 2.1(b), we have 

a(a;o)^^(a:o) 


and then 


u{xo)(^9 - u{xo) - Y 
u{xo) < 9 — 


+ mu{xo) + kv{xo) 
a{xo)v{xo) 


> 0 , 


1 + mu{xo) + kv{xo) 


< 9 . 
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Due to Lemma 2.2, we arrive at 0 < u < 0 on D. Similarly, we can show /r_|_ < v < 
+ l+mf+k^+ ^1- 

The (7^’“ boundedness of solutions comes from the elliptic regularity theory together 
with the Sobolev embedding theorem. □ 

We will use the local bifurcation theorem of Crandall and Rabinowitz [3] and the global 
bifurcation theorem of Rabinowitz [24] to prove Theorem 2. 

Denote the semitrivial solution curves by 

r„ = {{e, u, v) = { 0 ,o,^^)■ 9 > 0}, r, = {{0, u, v) = {0,0, 0); 0 > 0}. 

Define 


X = w^’P{n) X w^’P{ni), Y = LP{n) x lp(Di) with p>N, z = cl{n) x c^(Di), 


where 

= {w€ w^’P{ny, ^ = 0 on 5D}, cim = {w£ c^(ny, ^ = 0 on 5D}. 

The Sobolev embedding theorem implies X Y Z. 

Let solve 

A(j)* + {9* - = 0 in D, ^ = 0 on 9D, 

Alp* — pLip* -\ -— (j)* = 0 in Di, = 0 on dfii. 

1 + kp on 

Then^* = (-A + M/)S;Tffe'#'‘- 

Proposition 2.1 Let fj, > 0. Then there are positive solutions of (1.3) bifurcating from T^ 
if and only if 9 > 0 *{p,Qq), possessing the form 


Ti = {{0,u,v) = (6l(s),S(/>* + o(|s|),/r + sV’* + o(|s|)); s e (0,cr)} (2.2) 

with (0(0), u(0), v(0)) = (0*,O,/i) for some a > 0 in a neighborhood of (9*,0, p) G M x X. 
Proof. Denote by R = v — p, 


F(0,u,V)=(^ 


Au + fi(9,u,V + fi) 

AR + f2(ti,u,V + p) 


and Fi(9, u, v) = 


Au + fi(0,u,v) \P' 
Au + f2(fJ.,u,v) 


(2.3) 


with 


fi(9,u,v) = u(0 -u- 


a(x)v 


-), f 2 (h,'u,v) = v(p-v + 


cu 


1 + mu + kv 1 + mu + kv 

Obviously, F(0,u,V) = 0 is equivalent to Fi(0,u,v) = 0, and Fi(0,O,pi) = F(0,O,O) = 0 for 
0 G M. A direct calculation yields 

rKy)(«.0.0)[^.,^l = ( 

V Aif - pif + 


(2.4) 
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By the Krein-Rutman theorem, y)(0,0,0)[(/),-0] = (0,0) has a solution (/> > 0 if and only 
\i 6 = 9*. So (0*,O,//) is the only possible bifurcation point from which positive solutions of 
(1.3) bifurcate from r„. Besides, we have 


KerF(„y)(6l*,0,0) = Span {((/>*, -0*)}, dimKer F(„y)(r, 0, 0) = 1. 


For (</>, V’) G Y n Range F(„ y)(0*, 0, 0), choose ^ ^ such that 


A0 + (e-^W = 0. 


Aijj — + 


1 + /c/i' 
CIX 

1 + /c/r ' 


= Ip. 


(2.5) 


Multiplying by p* on both sides of the first equation of (2.5) and integrating by parts over 
n, we get pp*dx = 0. Then 


Range F(„y)(6»*, 0,0) = G Y; J^pp*dx = o|, (2.6) 

and thus 

codim Range F(„ y)(0*, 0,0) = 1. 

By a simple calculation, 


F0(r,O,O) = Fee(r,0,0) = (0,0), 

Fe{u,v){G*,0,0)[p*,p*] = {p*,0) i Range F(„y)(6i*,0,0). 

In conclusion, the proposition is proved by the local bifurcation theorem [3]. □ 

Proposition 2.2 Let — ^ < /r < 0. Then there are positive solutions of (1.3) bifurcating 
from r.u if and only if 9 > 9^. = having the form 

Fs = {{9,u,v) = {9{s),9 + sp^{x) + o(|s|),s + o(|s|));s G (0,d-)} (2.7) 

with 9{0) = p* = {A — g/)-i for some a > 0 near {9,9,0) G M x X. 

Proof. Let w = u — 9, 


Au; + (u; + 9){-w - x t 

Lr[d,w,v) - [ A , ^ , c{w+e) ^ j > 

V Av + v{fi - V + J 


Then Fi{9,u,v) = 0 is equivalent to G{9,w,v) = 0. We have 

G{w,v){G,w,v)[p,p] = 

Ap — {2w + 


a{x)v 


l+m{w-\-0)-\-kv 
c(w^6) 


' + 


a{x)mv{w^0) 
(l+m(Ri+0)+/cp)^ ' 

ckv{w-\-9) 


a(x){w+e){l+m{w+e)) j 
(l+m(L(;+0)+fci;)^ ^ 

ci;(H-fcp) 


+(^ - 2u) p - ^^ 


( 2 . 8 ) 
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T 


Ge{0,w,v) 


a{x)v . a{x)mv{w-\-0) 

l+m(if;+0)+/ci? ' {l-\-m{w^6)-\-kv)'^ 

cv{l-\-kv) 

{l-\-m(w+9)-\-kv)'‘^ 


The equation 0,0)[(^, ■0] = (0,0) is equivalent to 


A(l)-d 4 >- = 0 in O, 

< A^p + fiip + J^^p = 0 in Oi, 

= 0 on dQ, = 0 on (90i. 


(2.9) 


The second equation of (2.9) has a solution > 0 if and only if // = — , i-e. 0 = 

— c+m/^ ~ Thus { 6 ^, 6 , 0 ) is the only possible bifurcation point along T^, and solves the 
first equation of (2.9) with 6 = 6^ and ip = 1. It is easy to verify that 


KerG(^„^^)(6'*,0,0) = Span{(()>*, 1)}, dim KerG(^„^^) (0,0,0) = 1. 


A direct calculation shows 


G0(0*,O,O) = G00(0*,O,O) = (0,0), 

Range G(^„^^)(0*, 0,0) = {{f,g) €Y; gdx = 0}, codim Range G(^„^^)(0*, 0,0) 

Jn 


= h 


a(x) 


^6»(«),i;)(^*)0,0)[(/>*, 1] ( 0* ^ ^ ^ RangeG(^yj yj(6^,0,0). 

By the local bifurcation theorem [3], we get the desired results of the Proposition 2.2. □ 

In order to use the global bifurcation theorem for g > 0, define F 2 : M. x Z —>■ Z by 


^ ’ v-g ) V (-A +- g +f 2 {g,u,v)) 


( 2 . 10 ) 


Then (1.3) is equivalent to F 2 { 6 ,u,v) = 0. Let Ti C M x Z be the maximal connected set 
satisfying 


Ti C fi C {(0, u,v)eMx Z\{{e*,0, g)}-, ^ 2 ( 0 , u, v) = (0,0)}. 

From the global bifurcation theory of Rabinowitz [24] , one of the following non-excluding 
results must be true (see Theorem 6.4.3 in [18]): 

(a) Ti is unbounded in M x Z. 

(b) There exists a constant 9^9* such that {9,0, g) G fi. 

(c) There exists (0, ^,'tp) (Yi\{(0, g)}) with Yi = {{(p, Pp) G Z; (p(p*dx = 0} such that 

{9,(f),Pp) G fi. 


9 












Now we give the proofs of Theorems 2 and 3. 

Proof of Theorem 2. At first we know that u,v > 0 for any (9, u,v) € Ti which means that 

the case (c) above cannot occur by 4>* > 0. Otherwise there is a {6,u,v) G Ti such that (1) 

u > 0 with ^(xo) = 0 for some xq G Oi, or (2) u{xi) = v{x 2 ) = 0 for some xi G O and X 2 G Oi, 

or (3) h > 0 with u{x^) = 0 for some xs G VL. Denote by = {(/> G C^(D); 0 > 0 on D}. 

Choose a sequence {(^i, Uj, C Ti n (M x x such that lim (0^, Ui,Vi) = {6, u, v) 

_ 2—)-00 _ 

in M X Z, where 6 can be oo. Obviously, {u,v) is a non-negative solution of (1.3) with 6 = 6. 
By Lemma 2.2, one of the following must hold: 


(1) tt > 0, u = 0; (2) u = 0,v = 0; (3) u = 0,v > 0. 


For (3), we have —Av = v{fi — u) in Di, = 0 on 9Di, and thus v = ^. By Proposition 
2.1, this implies 6 = 6*, a contradiction to the definition of Fi. 

Suppose (1) or (2) is true. Integrate the second equation of (1.3) on Di with {u,v) = 
{ui,Vi) to obtain 



CUj 

1 -|- mui + kvi 


)dx = 0, 


i G N. 


On the other hand, /r > 0 and u = 0 ensure /r — Uj > 0, and thus fi — Vi+ > 0 for z 

large enough, also a contradiction. 

The case (b) is excluded by Proposition 2.1. So, the only true case is (a). 

From Lemma 2.3, {u,v) are uniformly bounded in Z as {6,u,v) G Fi which shows that 
6 is unbounded. Combining this with Proposition 2.1, we know that (1.3) has at least one 
positive solution for 6 > 0*(/i, Do) with fi > 0. 

Now, let {u,v) be a positive solution of (1.3) with m < • A direct calculation 

yields that u -|- By the monotonicity of the eigenvalue, we conclude that 


0 = Ai( — 


aix)v X ^ ^ (-0 + 

1 + mu + kv ^ 1 + kfi' 


Then 

0>Ai(:^^) = r(/r,Do). 

1 + kfi 

This shows that (1.3) has no positive solution whenever 6 < 6*{fi,^lo) and m < 
□ 


Proof of Theorem 3. When /r = 0, fix 0 > 0. By Lemma 1.1 and Theorem 2, we can take 
a sequence {(^i, u*, Uj)}“^ such that {ui,Vi) is a positive solution of (1.3) with fi = fii > 0, 
limj_>.oo fii = 0. By Lemma 2.3 and embedding theorem, we can choose a subsequence (still 
denoted by {{fii, Ui, u*)}^]^) such that {ui, Vi) converges to (ft, v) G Z, a non-negative solution 
of (1.3). By Lemma 2.2, u > 0 or u = 0 in D; u > 0 or {; = 0 in Di. 

If u = 0 and i; > 0, then fii — Vi+ , < 0 in Di for i large enough. This contradicts 

fn, - D + T+^^^Jdx = 0. 
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If it > 0 and D = 0, then — Vi + > 0 in fli for i large enough, also a 

contradiction with Viim - Vi + = 0. 

If ix = 0 and 1; = 0, then 6 — Ui+ > 0 in fl for i large enough, a contradiction 

to u^{9 -Ui + = 0. 

In summary, we must have u, u > 0 in H and Hi, respectively. This means that (1.3) 
possesses positive solutions for all 0 > 0 if /X = 0. 

Now, suppose —^<fi<0. For 6 > — > 0, the existence of positive solutions can 

be obtained from Proposition 2.2 by a similar global bifurcation analysis of F^ as that for 
the branch F„ with /x > 0. We omit the details. 

Conversely, let {u,v) be a positive solution of (1.3) with /x € (—^jO]. Then 0 < u < 9 
by Lemma 2.3, and hence 


/X = Ai(u - 


cu 


1 + mu + kv 


) > Ai(- 


cu 


1 + mu 


)>Ai(- 


c9 


1 + m9' 


c9 

1 + m9 ’ 


namely, 9 > 


□ 


3 Uniqueness of positive solutions 

In this section, we use topological degree to prove Theorem 4 for 0 > 9i and large /x. At 
first, introduce an auxiliary problem 


—Axx = ul 9 — u — 


a{x)v 


—Av = — V + t 

du 0 ^ 

dn dn ’ dn dfii 


1 + mu + kv 
cu 


1 + mu + kv 
= 0 


in n, 
in All, 


(3.1) 


with the parameter t G [0,1]. Eq. (3.1) reverts back to (1.3) if t = 1. When t = 0, we have 
from the second equation of (3.1) that v = fj,, and then obtain the scalar problem 


—Axx = u{9 — u — 

du 

dn 


a(x)/x 
1 + mu + 



in 12, 
on 912, 


which yields Eq. (1.6) as /x —>■ oo. 


(3.2) 


Lemma 3.1 Problem (1.6) has a unique positive solution if and only if 9 > 9i. 


Proof. Suppose 9 > 9i. Let (/> > 0 be the normalized eigenfunction with respect to 9i. Set 
XX = ecj). Then 


-Au = -eA(j) = e{9i - go(a^))</> = - qo{x) - ecp) + e0(0i - 9 + ecf). 
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Choose e small enough such that 9i — 9 + ecj) < Q to get 


du 

—Au < u{9 — qo(x) — u) in Q, = 0 on d^l. 

on 


Obviously, u = ecj) and u = 9 are a pair of positive sub- and supersolutions of (1.6) with 
u < u. We can get a positive solution of Eq. (1.6) by the sub-supersolution method. Let u 
and u be the minimal and maximal positive solutions to (1.6), respectively. Since 

/ Vu • Vudx = / uu{9 — u — qo{x))dx = / uu{9 — u — qQ{x))dx, 

Jn Jn Jn 

/ uu{u — u)dx = 0. 

Jn 


we conclude 


Therefore u = u. 

On the other hand, it is obviously true for any positive solution ui of (1.6) that 9 = 
+ qo{x)) > Xi{qo{x)) = 9i. □ 

Next, we show the uniqueness of positive solutions to (3.2). 

Proposition 3.1 Suppose 9 > 9i. There is a jl = fl{9) > 0 such that for any p, > ft, problem 
(3.2) has an unique positive solution. 

Proof. Since —q(){x) < — +kfi ’ then Ug^g^ is a subsolution of (3.2). Obviously, 0 is a 

supersolution of (3.2) and Ug^qg < 9. Then there exist positive solutions to (3.2). 

To prove the uniqueness of the positive solutions to (3.2), we at first show that the 
positive solutions of (3.2) are linearly stable for large p. Let U he a positive solution of (3.2). 
Consider the eigenvalue problem 


-A(t) = 9(t)- 2U(p - = 0 on dfl 


(1 -I- mil + kp) 

with the principal eigenvalues denoted by 


T] = r](p) = inf [ [|V())p — 9(j)^ + 2U((J + 


dn 




(3.3) 


(1 -|- mil + kpY 


(3.4) 


We have 


O = Ai(-0 + 2C/ + y^^^ML±^_^) > Xi{-9 - p) =-9 - p, 


(1 -I- mU + kp) 

i.e., rj > —9. Denote by p* the principal eigenvalue of the problem 


— A(j) = 9(1) — 2Ug nn{x)(l) — qoMf) + p*4> in D, = 0 on dVl 

on 
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(3.5) 






with the normalized eigenfunction 4>* > 0. Then t/* = jp ^ i+mu+kfi 

qoix) uniformly on fl as /r —>■ oo, we know U —)■ Ue^qg^x) uniformly on fl. It follows from (3.4) 
that 


V < 


6*^ + 2170*^ + 


a{x)n{l + kn) 2i, 
j(p \dx 


,• +1 mu - Us,,) + 


(3.6) 


Thus —6 < r] < M with M > 0 independent of fj,. We claim that liminf^_^(X) r] = r > 0. In 
fact, choose a sequence —)■ oo such that r]n ^ r, and 

AA -QA r, , a{x)nnil + knn) , , , ■ o _ n On /m 

^Yn — (^Yn ^'^nYn /-i , ,7 \nYn 'HnYri — 0 v^*') 

(1 + mUn + KAin) on 

with normalized (pn > 0, i.e. ||</>n ||2 = 1- As \Vpn\‘^dx are uniformly bounded with respect 
to n, there exists a subsequence cpnk 4>o weakly in H^{Q). Obviously, (/)o > 0 and ||<^o ||2 = 1- 
Multiply (3.7) by G Cpp{Q) and integrate by parts to have 


Vpn ■ 'Vifdx = 


[ [Gpr 

Jn 


(p 2Un,(pnP 


a(a;)/in(l + kpn) 

(1 + mUn + kpriY 


PnP + rjnPnPldx. 


Since Un Ug^qq uniformly on O as n —>■ oo, we have 


V(^0 ■ V<y9(ix = 


/ [d<Poq^ 

Jn 


- 2Ug^qqpop - qo{x)(j)op + rpo^jdx. 


Comparing with (3.5), we prove the claim that r = rj* > 0. So, there exists fl > 0 such that 
r] = r]{p) > 0 when p > jl, which implies the linear stability of the positive solutions to (3.2). 
Let 

HA u) = [M7 - A1->(M + « - „ - 
A = {tt G (7(11); eq < u < 6 + 1} 


with 0 < eo < xaiii^^^Ug^qq{x), 0 < t < 1, M large. Dehne 

S{t, u) = u — H(t, u). 


It is easy to see that 5(t, u) for all u G dA, 0 < t < 1. For large M, by the compactness 
of H, there are only finitely many isolated fixed points in A, denoted by rti,..., Um- Together 
with the linear stability of the positive solutions and the homotopy invariance of fixed point 
index, we have 

m 

1 = index(5(0, u). A, 0) = index(5(l, u). A, 0) = index(i7, Uj) = m. 

i=l 
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Therefore, there is an unique positive fixed point to (3.8) with t = 1 whenever fi > jl, i.e. 
problem (3.2) has an unique positive solution. □ 

Now, we can deal with the uniqueness Theorem 4. 


Proof of Theorem 4. Let {u,v) be a positive solution of (1.3) with large fi. Linearize the 
eigenvalue problem of (1.3) at (u, u) to have 


—Acp = 6(j) — 2u(p — 
—A'tp = — 2 u'i/; + 


a{x)v{l + kv) a{x)u{l + mu) 

(1 + mu + /cu )2 
cv{l + kv) 


(1 + mu + kv)"^ 
cu{l + mu) 


?/) + 


d(t> n n 

— = 0 on oSZ, = 0 on aili. 

on on 

Here (j), '0 and rj may be complex-valued. 

From Kato’s inequality, we have 


(1 -|- mu + kv)^ (1 + mu + kv)"^ 


'll) + r]4> in H, 
(p + rfip in Hi, 


(3.9) 


-A\(j)\ < -Re 



= Re( 0 | 0 | — 2u\4>\ — 
< 9\(j)\ — 2u\(j)\ — 


a{x)v{l + kv) 
(1 -|- mu + kv)"^ 
a(a:)u(l + kv) 


a(x)u(l + mu) , (p , 


(1 -|- mu + fcu)^ 


+ 


a{x)u{l + mu) 
(1 + mu + kvY^ 


+ Re(r/)|<).|. (3.10) 


To obtain the linear stability, it suffices to prove that for any h > 0, there exists > 0 such 
that the eigenvalues t] of (3.9) satisfy Re( 7 /) > rj* — 6 when > ns- Otherwise, there exist a 
ho > 0 and a sequence {(/i„,r/„,</>„, V’n)}^! satisfying (3.9) with \\(pn\\2 + WYnh = 1) 
and fin ^ oo as n ^ oo such that Re(r 7 „) < rj* — ho- Replace {fi,r],u,v,(pA) ia (3.10) with 
{fJ‘n,iln,Un,Vn, (pnAn), multiply by \(pn\, and then integrate by parts over H to have 


/ 

Jn 


\V\4>n\\‘^dx < [ {e\(pn\‘^ - 2Un\(pn\‘^ - 
Jn 

a{x)un{^ + mu 


a{x)Vn{l + kVn) 
(1 -I- mUn + kVn) 




. , , . 2 \Yn\\(pn\)dx + {if - 6o) \(pn\‘^dx. (3.11) 

(^i mUn ~r HVn) Jn 


Let Vn be the principal eigenvalue of the eigenvalue problem 

A /) o a{x)vn{l + kVn) , ■ n n ac 

-A(p = Bip- 2un^p - 7 :—-+ rn(p m H, = 0 on Oil. 

(1 + mUn + kVn) on 


We know that 


He 

rn-rf 


mi - 


9 A + “JUnA + 


a{x)vn(l+kvn) 2 
{l+mUn+kVny-‘ ^ 


In ‘^'^dx 


r]*A]dx 


and Vn —)■ rj* by the proof of Proposition 3.1. So, there exists a W > 0 such that Vn — rj* > 
for n > N. Thus by (3.11), 


2 
















2 aij. |2 , o I , |2 , a{x)Vnil + kVn) ,, ,2 


— / [|^</'nl ~ ^|0n| +2ltn|0n| + 
JQ 


< —<5o 


[ \4>n\‘^dx+ f 
Jo. Ju 


(1 + mUn + kVnY 

aun{l + mUn) I , II , I , 

:\ll}n\\(Pn\dx. 


\^n\ - ri*Hn\ ]dx 


(1 + mUn + kVn)' 


Since ^ ™ as n ^ oo, then \cl)n\‘^dx 0 . 

Using Kato’s inequality again, we have 

Al/l^ 1/1 O l/li CUn(l + mUn) I ; I I + kVn) IJ, I i / * X W I I 

~^|wn| ^ /^nlwnl “ 2 Un,|wn| + T:j ' ' j w lU^I “ 1 “ T\ i To l^rnl \J1 ~ Oojlwnl- 

(1 + mUn + kVn) (1 + mUn + kVn) 

Multiply by |'0ri| and integrate by parts over fii to get 


/ W\'4^n\?dx< / [/inIV’nP - 2Un|V’nP + 

J Hi J ni 


CUnjl + mUn) 

(1 + mUn + kVnY 


{'(Pr, 


CVnil + kVn) I , II / I , / * CM/ |2i , 
+ 7W--—;- —\(Pn\\Vn\ + [1] -OojlUnl \dx 


< 


(1 + mUn + kVnY 

CUn(l + mUr, 


I [ IJ'n + 
J Hi 

/ 

JVi 


Oi (1 + fnUn + kVnY 

CVn{l + kVn) 


+ rf - doWipn^dx 


Consequently, 


[ I'lpnl'^dx < — [ 
jQl t^n Jn 


(1 + mUn + kVnY 

CUn{l + mUn) 


\4>n\\'>Pn\dx. 


1 

+ - 


(1 + mUn + kVnY 

CVn{l + kVn) 


< 


IJ-n (1 + rnun + kVnY 
I r , c9 


+ ri* - 6Q\\iin\ dx 


\(pn\\lpn\dx 


[ - 6o)\'lpnfdx + — / ^\(l)n\\P^n\dx. 

Jn, l + mO Hn Jn, k 


in 


l^n Joi 1 + m 6 » 

This concludes j'i/’npda; ^ 0 as n —>■ oo, since fin ^ oo and l^riMV^nl are bounded 
L 2 (Oi). 

In summary, we have obtained \4>n\'^dx,J^^ \'ipn\‘^dx —?> 0, as n —>■ oo which contradict 
with Wcpnh + llV'nIb = 1 - 

By using a similar argument to that in the proof of Proposition 3.1, we get from the 
linear stability of the positive solutions to (1.3) and Proposition 3.1 that the solution of (1.3) 
must be unique when fi > max{fl, fio} with fiQ = mi{fis;6 € ( 0 , 7 /*)}. 

Finally, we consider the asymptotic behavior of the unique positive solution (u, v) as 
// —)• oo. Since —)• 0 as /x —>• oo, for any e > 0, there is a /Xe > 0 such that 

i+rZ+kv < e for /X > fi,. Then 

fiv — < — Au < (/X + e)v — in fli, 

which yields /x < u < /x + e for /x > And thus u — /x—^-Oas/x^oo. We know that 
i+ml+kv 9o{x), and then u Ue^g^ix) uniformly on H as ^ oo. □ 
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4 Discussion 


In a reaction-diffusion system of predator-prey PDE model, in addition to the interaction 
mechanism between the species, the behavior of the species is also affected by the diffusion 
of the species, as well as the size and geometry of the habitat. Obviously, the prey species 
would die out under excessive predation from nature or humans. The results obtained in this 
paper show the way in which the created protection zone saves the endangered prey species in 
the diffusive predator-prey model with Beddington-DeAngelis type functional response and 
non-flux boundary conditions. 

Compared with previous results on protection zone problems with other functional re¬ 
sponses such as the Lotka-Voltera type competition system [7], Holling II type predator-prey 
system [4, 9], and Leslie type predator-prey system [8], richer dynamic properties have been 
observed for the model (1.1) with B-D type functional response in this paper. It can be found 
that a total of four threshold values are obtained here for the prey birth rate 6, i.e., 9q, 9*, 
9i (for the predator growth rate n > 0) and 0* (for fj. <0). 

The first threshold value 9q gives the necessary condition for establishing a protection 
zone to save the prey u. By Theorem 1, the survival of u could be automatically ensured 
without protection zones whenever 9 > 9q = |, which can be realized when the refuge ability 
of the prey is properly large that k > or the predation rate is small that a < 9k. In other 
words, the protection zones have to be made only if the prey’s refuge ability is too weak 
with respect to its birth rate 9 and the predation rate a. This matches with the mechanistic 
derivation for the B-D type functional response proposed in [12] . In addition. Theorem 1 says 
also that the protection zones are unnecessary if the predator’s growth rate /i < 0, where the 
predator species v can not live without the prey u, and thus the extinction of v cannot take 
place after of u. 

The second threshold value is 9* = 9*{fi, Dq) = Ai(g(x)) with q{x) = and /r > 0. By 
Theorem 2, the positive steady states can be attained for 9 € {9*, 9q). Due to the monotonicity 
of the principal eigenvalue Ai = Ai(g(x)) with respect to q{x), we know that the threshold 
value 9* would be enlarged (and hence harmful for the prey u) when the predation rate 
a{x) or the predator’s growth rate /i increase, or when the prey refuge k or the size of the 
protection zone Dq decrease. Conversely, Theorem 2 says also that the prey u must become 
extinct when 9 < 9* with the handling time m of vh being shorter than All of 

these match with those in [12]. In addition, since 0*(/x,Do) < 9q is strictly increasing with 
respect to fi and decreasing when enlarging Dq, letting /r —>■ oo, we get the third threshold 
value 9i = 0i(Do) such that if 9 > 9i{Qo), prey species could survive no matter how large the 
predator’s growth rate is. The critical 9 = 0i(Do) implies a critical size of the protection zone 
as well, namely, if the real protection zone Dq D Dq, the survival of the prey with such birth 
rate 9 is independent of the predator’s growth rate. Also, the uniqueness and linear stability 
obtained in Theorem 4 for ^ large enough are reasonable because —)■ 0, and hence 

u — //—>-0as/r^oo. 

Since the condition fi <0 yields the survival of u without protection zones by Theorem 
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1, the fourth threshold value 0* obtained in Theorem 3 with // € only made for v 

alive. In fact, the conversion of prey is limited by as shown in the proof of Theorem 1, 
the predator v must be die out if its growth rate // < — ^. With such non-positive growth 
rate fj, € (“m’ should be properly large number of prey to survive the predator, just 

as described via the criterion 0 > 0* = > 0 in Theorem 3. It is worth pointing out 

that the threshold value 0* for alive predator v would be enlarged as m (the handling time 
of vh) is increasing. This well matches the mechanism in [12]. 

We have shown the effect of refuge ability of the prey and protection zones on the the 
coexistence and stability of predator-prey species in the paper. In fact, protection zone 
can be regarded as another refuge offered by human intervention, which is necessary if the 
prey’s refuge ability is too weak in the predator-prey system to prevent the extinction of 
the prey populations. The critical sizes of protection zones, obtained in this paper and 
represented by the principal eigenvalues 0* = Xi{q{x)) and 9i = Xi{qo{x)), show the basic 
requirement (depending on the predator’s growth rate) and the sufficient one (working under 
any predator’s growth rate), respectively. The results of the present paper would be helpful 
to design nature reserves and no-fishing zones, etc. 
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